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Abstract 

The electronic properties of semiconductor, vertical, double quantum dot 
systems with few electrons are investigated by means of analytic, configuration 
interaction, and mean field methods. The combined effect of a high magnetic 
field, electrostatic confinement, and inter-dot coupling, induces a new class 
of few-electron ground states absent in single quantum dots. In particular, 
the role played by the isospin (or quantum dot index) in determining the 
appearance of new ground states is analyzed and compared with the role 
played by the standard spin. 

I. INTRODUCTION 

The behavior of a small number of electrons confined into a single quasi-two-dimensional 
quantum dot (QD) in the presence of a magnetic field has been studied over the past few 
years. 'Ehi Most of the work has focused on the high magnetic field regime of fully spin- 
polarized electrons where incompressible electron states analogous to the Laughlin states 
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of the fractional quantum Hall effect exist. Only very recently the problem of spin and its 
implications has been addressed. ih§ The ground and excited states of a spin-unpolarized 
QD turn out to be much more complex than those of the polarized QD. The role of spin 
and spin induced interactions can be effectively simulated by the isospin (layer index) in 
double layer systems where the isospin-up isospin-down states correspond to electrons on 
layer 1 or layer 2. Recent experimental iHll and theoretical studies of a double 

layer two-dimensional electron system (DL2DES) in a high magnetic field have shown a 
rich variety of new incompressible states related to the quantum Hall effect il (QHE). For 
instance, in recent experiment of Ref. 23 S. Q. Murphy et al. have reported the existence 
of a new incompressible QHE state at a filling factor of individual layers z/ = 1/2 §1 in the 
absence of inter-layer tunneling. This coupling-induced QHE state has no analogue in single 
layer systems. The distance between layers, i.e., the strength of the Coulomb coupling, 
determines whether the incompressible state appears or not. This experiment seems to 
confirm the predictions made in Refs. 15-17, and 21 that the new QHE states are supported 
by the inter-layer Coulomb interactions and by the single-particle symmetric-antisymmetric 
gap in the presence of arbitrary tunneling. The breakdown of these states occurs when the 
Coulomb coupling is switched off and the intra-layer correlations become dominant. El^il 

In this work we examine the role of confinement, magnetic field, and inter-dot coupling in 
the generation of a new type of few-electron states in a vertical double quantum dot system 
(DQDS). We concentrate on the high magnetic field limit where electrons in the ground 
and low-lying states are spin polarized due to Zeeman energy, and the effect of the isospin 
degree of freedom can be isolated, the appearance of several ground states (GS) with an 
unexpected isospin is the most striking feature of such systems in comparison with isolated 
polarized QD's . 

The paper is organized as follows: Section I is devoted to the description of the model; 
in Sec. II we present the analytical results for the simplest possible case of two electrons 
in a DQDS; in Sec. HI we present numerical results for up to six electrons, and Section IV 
contains conclusions of the work. 



II. THE MODEL 



We consider a pair of identical vertical coupled QD's with parabolic in-plane confinement, 
containing electrons, and separated by a distance D as shown in Fig. |l|. The normalized 
single-particle states \m,n\a > of each individual QD (in the presence of a magnetic field 
B oriented normal to the plane of the QD) are simply harmonic oscillator states 

\m, n- a >= ./^(at)-(6t)"|o, 0; a > (1) 
V mini 

with a labeling the QD index or z-component of the isospin quantum number. In analogy 
with the standard spin, a takes on the values +1/2, —1/2 (for brevity, +, — , from now on) 
for electrons in the "upper" and "lower" QD, respectively. The single-particle energies are 
those of a pair of harmonic oscillators (the Zeeman energy is omitted and h = 1 for the rest 
of this work) 

e^na = ^+{n+^) + ^-{m+^) (2) 

with = [\J ujI + Aujq ± a;c]/2. uJc is the cyclotron frequency and ujq is the frequency 
characterizing the parabolic confinement of both QD's. We shall restrict to the case of high 
magnetic fields (5 > 2 — 3 Tesla) where electrons are spin polarized, and the condition 
ujc > ujq is satisfied. In this regime one can restrict the basis set to n = and the problem 
becomes one-dimensional. We shall omit index n in what follows. 

Allowing for inter-dot tunneling, the full-interacting Hamiltonian of the DQDS can be 
expressed in second quantization as 

H = J2 \^<y.<y'^maC^ma^ma' + (1 " ^ a ,<y')iC^ma(^mA + 
man' 

1 \ " T/mim2m3m4 t ^ p ,p 

where c"*", c are the single-particle creation and annihilation operators, \/j^i™2m3m4 
Coulomb interaction terms, and t is the inter-dot hopping matrix element. We have in- 
cluded here inter-dot direct and exchange terms but neglected some off-diagonal scattering 
elements, negligible in the weak hopping (t 0) limit. 
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Alternatively, we can use an isospin 1/2 operator representation to define isospin opera- 
tors in terms of previous creation and annihilation single-particle operators as 
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Omitting linear terms in Pmim2 and Cm2m3 allows us to emphasize the appearance of new 
isospin-isospin interactions in the DLQD Hamiltonian: 

rr _ V^/ \f/-x \ \\ \^ T/mim2m3m4 

m mim2m3m4 

I ^ T^mim2m3m4 >z /a\ 

^ m,\m2ra'ira4 

The Hamiltonian (|) contains two S\J (2) symmetry-breaking fields: The first one along the 
X direction and proportional to the hopping t, and the second one along the z direction, 
the isospin-isospin interaction, proportional to Ve (in addition to the usual charge-charge 
interaction). The Coulomb matrix elements are given by 

V'mim2m3m4 i-\Tm\m2'm,-im4 T/mim2m3m4\ 

D — 2^ ++ +- ' 

V'mim2m3m4 /T/mim2m3m4 rm\ni2mzm4\ 

E — 2^ ++ ~ +- 

where (identical to V ) are intra-dot Coulomb matrix elements, and V^+_ are inter-dot 

Coulomb matrix elements. 

As far as the symmetry-breaking term t is concerned, two physical situations can be dis- 
tinguished: the "incoherent" and the "coherent" one. The incoherent case {t = 0) describes 
isolated QD's coupled only by Coulomb interactions. Electrons cannot transfer (exchange) 
between QD's, and are localized on individual QD's. This distinguishability manifests itself 
in the anticommutation relation {c^^g., 0^20-'} = Sa,a'Smi,m2- This is the well known "layered 
electron gas" model. 

The second case is the coherent behavior of QD's where electrons cannot be in a specific 
QD but occupy the symmetric (s) or antisymmetric (as) orbitals of a pair of QD's. The 



transformation from orbitals localized on individual dots to the symmetric/antisymmetric 
orbitals is equivalent to the rotation of the isospin. Let us define a rotated isospin repre- 
sentation, {a}, which diagonalizes the hopping part of the Hamiltonian (^) for arbitrary 
hopping matrix element t: 

(5) 

We wish to emphasize that this transformation does not depend on the strength of the 
hopping matrix element. We can define isospin operators p, in the space of coherent 

operators {a} to write the coherent Hamiltonian as 

m mim2m3m4 

I ^ \ " T^mim2m3?Tt4 >a; /r.\ 
' / J >>m\m4'^m.2va-i' y^l 

^ mim2m3m4 

In this rotated isospin space the hopping matrix element t is simply equivalent to an external 
field. The coherent Hamiltonian (|^) is similar to the previous one in the sense that it 
also presents two SU (2) symmetry-breaking terms: that proportional to the hopping t and 
the isospin-isospin interactions proportional to Ve (again, in additon to the normal SU{2) 
invariant charge-charge interactions). This isospin-isospin interaction Ve, present in both 
Hamiltonians, is, ultimately, responsible for the novel physics in a DQDS. Since this physics 
is dominated by Coulomb interaction we shall concentrate on the t — limit. 



III. ANALYTICAL RESULTS FOR THE TWO-ELECTRON DQDS 

The case of = 2 is the simplest case that deserves to be studied in detail and can be 
solved analytically. We can write the total wave function \E'(ri, zi r2, Z2) as a product of 
the center of mass (R) wave function, relative motion (r) wave function, and the rotated 
isospin wave functions (symmetric or antisymmetric) 
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*(ri,r2) = V'''^'"(i)0^'^(r)|i,j >= \M'^,M^;i,j > i,j^s,as. (7) 

The center of mass motion (characterized by an angular momentum M'^) separates from 
the Hamiltonian, and isospin states can be characterized by 4 orthogonal rotated isospin 
states {\s,s >, \as, as >, \s, as >, \as, s >}. These states can be written in a more familiar 
language of the usual isospin states, e.g., {|+, + >= |/ — l,Iz — +1 >}, based on orbitals 
localized on individual QD's: 

\s,s> = ^ { l^=^'^-=+^X^'^-=-^> + |7 ^ 1,4 ^ >} 

\as,as> = ;j5{^-i^-±i^^-i^-|/ = l,4 = 0>} 

\s, as> = ^ {|/ = 1, 7, = +1 > -|7 = 1, 1, = -1 >} 

\as,s > = |7 = 0, = > . 

The first three states correspond to a well-defined isospin I — 1 but an undefined z com- 
ponent of I. The expectation value of for the three coherent states 7 = 1 is zero, but 
quantum fluctuations are present in the 7 = 1 states, in contrast with the incoherent 7 = 
state. This means that the two 7 = 1 states -^{\s,s > +\as,as >} and \s,as > are not 
eigenstates of 7^, and one cannot determine on which QD the electrons are locahzed. They 
correspond to having both electrons in one or the other QD of the DQDS. On the other 
hand, the |7 = 1, 7^ = > and |7 = 0, 7^ = > states (i^lls, s > —\as, as >} and \as, s >, 
respectively) correspond to electrons on opposite QD's and their energy is determined only 
by inter-dot interaction V+- = Vd — Ve- 

We can now write the relative particle Hamiltonian in the rotated isospin space as a 4 x 4 
matrix. It is easy to see that the first three states correspond to total isospin 7=1 and 
hence the in-plane relative particle wave function must be antisymmetric. This corresponds 
to odd relative angular momentum, M*", of the relative particle. The Hamiltonian 77i for 
the 7 = 1 states can be written as 
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t + VL_M-'+ < M^\Vd\M^ > <M''\Ve\M^> 

<M''\Ve\M^> -1 + ^1^^+ < M''\Vd\M'' > 










and for / = we simply have 

n_M''+ < M''\Vd-Ve\M'' >= Ho, (10) 

with relative angular momentum M*" being even. 

From the rotated isospin Hamiltonian it is clear that (a) only the symmetric- 
antisymmetric states are coupled, (b) coupling is due to the symmetry-breaking exchange 
interaction Ve, (c) the coupling between symmetric- antisymmetric states is present even in 
the absence of tunneling (t = 0), and (d) only these coupled symmetric-antisymmetric states 
are affected by inter-dot hopping (t). 

The two-electron Hamiltonian can be easily diagonalized and simple analytical expres- 
sions for energies and wave functions obtained (not shown here for brevity). Figure ^ shows 
the evolution of the total energy spectrum as a function of the distance between QD's for a 
given value of the total angular momentum M = ikf + M*^™- = 6 and zero hopping matrix 
element t. Standard values of dielectric constant and effective mass for GaAs have been 
taken throughout the calculations and the confining energy of the QD's in the DQDS is 
taken to be 5 meV. From now on it is also convinient to define an inter-dot coupling con- 
stant: a = One can see in Fig. ^ the energy splitting of the 1 = 1 states (when 
odd), which are degenerate when a = 1, as the distance increases (i.e., as the coupling 
constant a lowers its value). Those with = go down in energy while those with = ±1 
remain degenerate and constant (notice that this would not be the case any more if t 7^ 0). 
It must be pointed out that the total isospin J is a good quantum number for any value 
of the coupling constant in the two-electron DQDS. This is no longer the case for a higher 
number of electrons as will be shown below. In Fig. ^ we show the phase diagram (isospin 
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/ and angular momentum M) of the two-electron DQDS GS as a function of the magnetic 
field B and coupling constant a. The magnetic field changes the ratio of kinetic {^-) to 
Coulomb energy and induces changes in the isospin and angular momentum of the DQDS 
GS. Similar transitions of the isospin and angular momentum of the GS can be induced by 
changing the coupling constant as shown in Fig. ^ 

IV. MORE THAN 2 ELECTRONS IN THE DQDS 

A. Zero distance limit: No SU(2) symmetry-breaking interactions 

The zero distance limit (ZDL) (D = 0, a = 1), i.e., the limit of having the QD's forming 
the DQDS superimposed in real space, although urealistic, is interesting and deserves to be 
studied in detail. In this particular case (always with t —>■ 0) the SU{2) isospin symmetry 
is not broken by the isospin-isospin interactions since Ve = 0. In the ZDL the DQDS of 
identical QD's under high magnetic fields is completely equivalent to a single QD with Lande 
factor g ^ 0, i.e., in the zero Zeeman limit. The role of spin in the DQDS (frozen out by 
the magnetic field) is now played by the isospin. In the case of non-identical QD's, the 
presence of an "isospin Zeeman energy" (for instance, a difference between the confinement 
energies of each QD) would make the DQDS equivalent to the spin-polarized single QD . 
The feasibility of fabricating identical QD's with zero "isospin Zeeman energy" is one of the 
most appealing possibilities presented by such systems. As will be shown below, it is the 
fundamental origin of the new electronic properties that appear in a DQDS at high magnetic 
fields, compared to those appearing in the same regime of fields in a single QD (i.e., in a 
spin-polarized QD). 

We now extend our study to a larger number of electrons by employing exact diagonaliza- 
tion techniques. The intra- and inter-dot electron-electron correlations have been taken fully 
into account by expanding the many-body wave function in terms of Slater determinants 
(or configurations) and diagonalizing the Hamiltonian (|^) for 4 and 6 electrons. The Hilbert 
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space has been restricted to the states of the lowest Landau level of each QD {n = 0). Due 
to the circular symmetry, the diagonalization can be done in separate subspaces of configu- 
rations with the same z-component of the total angular momentum, M. In the absence of 
hopping any subspace of given M can be split, in turn, into orthogonal subspaces of given 
Iz- In this way the size of the matrices to diagonalize becomes smaller and computation- 
ally more accesible. The GS's and lowest-lying eigenvalues and eigenfunctions of all the 
subspaces of different M were obtained using standard diagonalization routines. 

Figures ^(a) and (b) show the evolution of M for the DQDS absolute GS as a function 
of Q- (i.e., B) for 4 and 6 electrons . The absolute GS angular momentum M goes through 
a series of increasing values as lowers its value (i.e, the magnetic field rises). The 
competition between kinetic energy and Coulomb repulsive energy determines the value of 
M: the kinetic energy (due to the confinement and magnetic field) favors electrons in the 
center of the QD, the Coulomb repulsion tends to spread the charge. 

The results in Figs, ^(a) and (b) show a remarkable stability of the GS's with M = 
6,1 = 2 for 4 electrons, and M = 15, / = 3 for 6 electrons, against changes in Q^. These 
stable GS's are the only ones appearing with maximum total isospin J, and are five-fold 
degenerate (/^ = 2, 1, 0, —1, —2) and seven-fold degenerate (/^ = 3, 2, 1, 0, —1, —2, —3), for 4 
and 6 electrons, respectively. In particular, the = states are described exactly in this 
ZDL by the Jastrow-type correlated wave function [111] @ 

(^{zi,...,Zn/2,Wi,...,Wn/2)=A[ Yl i^i-Zj) n (U^t-U^j) n iZi-Wj)x 

l<i<j<N/2 l<i<j<N/2 l<i,j<N/2 

N/2 N/2 

exp(- 5: |^,|V4/2-^ 1^,174/2) X 

i=l i=l 
1+ >i ...|+ >N/2 \- >1 --I- >7V/2], (11) 

proposed by Halperin in the context of fractional QHE wave functions with spin degrees of 
freedom, and used later in the context of a DL2DES. The electron coordinates of the 
upper dot are given by Zi = Xi — iyi and those of the lower dot by Wi = x[ — iy[, and -|- and 
— denote the values of the QD indices (see previous section). The symbol A represents the 
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antisymmetrization operator. Alternatively, these states can be also expressed in terms of 
the single-particle occupation numbers of each QD, z/^ and v~ . For instance, for 6 electrons 
z/^'~ = 0.5 for m = 0, 1,2,3,4, and 5, and for the following m's. If we define a filling 
factor for the DQDS as p^Q^^ = (z/+ + z/")^ where the brackets denote average over the 
lowest occupied m's, the GS's M = 6 (A^ = 4) and M = 15 {N = 6,) clearly correspond 
to jj^^^^ = 1. These p^^^^ = 1 states can be considered, in turn, as precursors of the 
incompressible state at total z/ = 1 observed in a DL2DES. 

In Fig. D we show the six-electron excitation spectrum of a DQDS with a p^^^^ = 1 
degenerate GS of angular momentum M = 15 (marked with a star in the plot). In partic- 
ular, the fully isospin-polarized state within the degenerate subspace has been considered 
as the absolute GS in what follows. The excitation spectrum consists of two branches: the 
branch with M < 15 and that with M > 15. Let us first concentrate on the branch of the 
spectrum with angular momentum M < 15. By examining the wave function we have been 
able to identify those excitations corresponding to isospin-flip excitations (solid dots in the 
plot). Those with M close to M = 15 can be associated to isospin-flip spin-wave-like edge 
excitations. @ Those with M farther below M = 15 correspond to isospin-flip quasiparticle 
excitations (magnetoexcitons il) and consist in flipping the electron's isospin and moving 
it from the edge of the u^^^^ = 1 droplet to a reversed isospin single-particle state m closer 
to the center of the DQDS. The value of the total isospin of such excitations is the maximum 
possible value according to the spin flip and to the subspace in which it is found, i.e, the 
value of M of the excitation (all those shown with solid dots in Fig. ^ on the left branch, 
correspond to / = (A^ — l)/2 = 5/2). We can see that isospin-flip excitations are not, in 
general, the lowest energy states in subspaces of given M, but there appear a few other 
states of lower energy and minimum value of / within each subspace. 10 

We have also identified isospin-conserving quasihole-quasielectron pazr excitations (open 
dots in the plot) consisting of a quasihole in a single-particle state of m < 5 and an electron 
added to the edge of the droplet. Similarly to the spin- flip excitations, those with M far 
from M = 15 present the better-defined character of a single, localized hole in the droplet. 
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These excitations present the maximum possible value of / {N/2 = 3). In contrast with the 
isospin-flip excitations, there appear many excitations of lower energy and minimum isospin 
within each subspace of given M. 

From the excitation spectrum we can understand the evolution of the absolute GS an- 
gular momentum with The stability of the u^^^^ = 1 GS's is nothing but the direct 
consequence of the cusp-like structure or gap exhibited by the excitation spectrum (see Fig. 

and such state can be referred to as an incompressible GS. Up to the values of Q- stud- 
ied, no other stable GS's seem to appear in our calculations. In the case of identical QD's 
almost all the lowest states of different subspaces of M will become the absolute GS of the 
DQDS at a certain value of B, which is in clear contrast with the case of a single, isolated 
polarized QD. In this ZDL the GS takes on almost all possible values of M for the range 
of variation of but such values change very quickly with Q_ (except the p^^^^ = 1), 
presenting no stability. The relevance of these kind of new absolute GS's in a single QD 
at low magnetic fields (when the Zeeman energy cannot be considered infinite) has been 
stressed in Ref. 8. Their importance lies in suppressing, through the spectral function of 
the system, the single-electron tunneling rates, and therewith, in strongly modifying the 
transport properties of the system. One would expect that the presence of such states in 
the DQDS would give rise to similar effects on transport properties of these systems as long 
as they do not disappear for a realistic situation, i.e, for a certain distance between QD's. 
Taking into account realistic values of the distance is the topic of the following sections. 

As mentioned previously, a significant difference between the confining energies of the 
QD's forming the DQDS (or the presence of the g factor for a single QD at high B) tends to 
favor fully isospin-polarized (spin-polarized for a single QD) absolute GS's. QD's with fully 
polarized electrons have been extensively studied over the past few years. It is known that, 
for instance, for four and six electrons, when only the spin-polarized states are relevant to the 
GS properties, the value of M for the absolute GS is restricted to a series of specific numbers: 
6, 10, 14, for 4 electrons 01 and 15, 21, 25, for 6 electrons 0. These numbers are 
known in the literature as "magic" numbers, i'i This can be seen easily from Fig. ^: The 
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lowest states in each subspace of M will correspond now (in the presence of isospin Zeeman 
energy) to isospin-conserving quasihole (M > 15) or isospin-flip quasiparticle (M < 15) 
excitations. Only those with quasiholes (quasiparticle) close to the center of the droplet will 
become the GS of the DQDS as B changes [those with downward-pointing arrows in Fig. 
^. A simple description of the magic absolute GS's (for M > 15) in terms of "bosonic" 
operators acting upon the electronic v = 1 droplet has been presented elsewhere, il 



B. Short distances: Weak SU(2) symmetry-breaking interactions 

We now discussed how the new minimum-isospin GS's of a ZDL DQDS evolve for realistic 
distances between QD's. Figures Ka) and (b) show, for D = lOOA and D = 50 A, the 
evolution with Q_ of the M value of the absolute GS (from now on we will restrict to the 
6 electrons case). As was clearly shown in Fig. as the distance between QD's increases, 
or the symmetry-breaking interaction term Ve is stronger, the isospin multiplet degeneracy 
is removed and the GS will always have Iz = 0, i.e., equal number of electrons in each dot. 
The total isospin, /, is no longer a good quantum number for D 0, but one can still trace 
it back to its original value aX D = 0, and use it to label the states as long as the distance is 
not too large. As can be seen in Figs, ^(a) and (b), for short distances {D = 50 A, a ~ 0.8) 
many of the absolute minimum-isospin GS's in the ZDL survive. The M = 15 GS also 
remains stable. On increasing the distance {D = lOOA, a ~ 0.6) the situation changes 
dramatically: The symmetry-breaking interactions have made all the minimum-isospin GS's 
of the ZDL disappear. Instead, many new stable states appear (M = 6, 9, 12, 15, 18, ...), and 
the only stable GS in the ZDL, M = 15, has become less relevant. At the same time, the 
overlap of such GS (or the Jastrow-type wave function shown above) with those of the same 
M at distances different from zero decreases with D. Instead of showing this overlap we 
have chosen to show in Fig. |^ the single-particle occupation numbers (z/~) in the ZDL, 
for D = 50A, and for D = lOOA. One can see that the regular occupations forming the 
^DQDS _ g^ate in the ZDL melt as the distance between the QD's increases, and the value 
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of M = 15 becomes less relevant for the GS. These new stable states have their origin in the 
"superposition" of the stable ones for each QD with 3 electrons (M = 3,6,9,...). i'i The 
next section will clarify what we mean by such superposition. 

C. Large distance limit: Strong SU(2) symmetry-breaking interactions 

We have seen in the previous section that, as the distance D between dots increases the 
situation seems to change noticeably. In order for this large-D limit to be understood we 
have carried out a a self-consistent Hartree-Fock (HF) treatment of the inter-dot coupling, 
but conserving the intra-dot correlation. The procedure is the following. A GS solution of 
total angular momentum, M+, by means of an exact diagonalization of, for instance, the full 
interacting Hamiltonian of the upper QD, ^ is found as described in the previous section. 
Then, the single-particle energies of the lower dot are modified by 

where, in the calculation of the self-energy S+, Hartree-Fock-like diagrams have been 
used. If the inter-dot hopping is forbidden, then only a Hartree-type diagram is allowed 
and corresponds to taking into account the single-particle occupation numbers of the upper 
GS, together with the coupling term \/^-i_™-2m2m.i ^ Now, an exact diagonalization of the 
lower full interacting Hamiltonian, if^, with the corrected single-particle energies is done, 
a lower GS of M~ is found, and the occupations of this GS are used to modify, in turn, 
the single-particle energies of the upper dot. The process continues until convergence is 
achieved. 

Figure shows the evolution of the renormalized GS energy of our DQDS as a function 
of the magnetic field (or VtJ) for the case of 3 electrons in each dot (J^ = 0) and for four 
different distances. The GS's for D = 200A (a ^ 0.3) correspond to a superposition of 
M+ = 3, M- = 3 (which gives us M = M+ + M" = 6), M+ = 6, M" = 6 (M = 12), 
and M+ = 9, = 9 (M = 18). The possible values for M+ and M" (3,6,9,...) are the 
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corresponding magic angular momenta M mentioned above for the case of 3 electrons in a 
single QD. HJi As we bring the dots together, these states shift to lower fields, so one can 
think of the inter-dot coupling in terms of an " additional" magnetic field which adds to the 
bare one to give a stronger effective value of B. 

Exact energies are also shown in Fig. ^ At large distances {D = 200A) both procedures 
give the same energy for the GS. As the dots are brought together, the exact GS energy 
becomes smaller than the HF energy. This fact is more noticeably for the shortest distances 
where this deviation takes place at B ^ 7 T and B ^ A T for D = lOOA and D = 50A, 
respectively. This deviation points out the fact that inter-dot correlations have grown to 
play their role, a role that cannot be described in terms of a mean field theory. A phase 
diagram is shown in the inset: The inter-dot correlation grows to be relevant with the 
inverse of the distance (as one would expect), and with the strength of the magnetic field. 
Now one can understand the appearance of the u^^^^ = 1 from a different point of view: 
The new incompressible state at short distances (M = 15) has its origin in the spatial 
inter-dot correlations, and cannot be obtained as a simple superposition of two equal QD 
stable electronic configurations. The origin of such particular occupation numbers can be 
understood. One can form a GS with M = 15 out of two single QD configurations with 
= 3 and = 12. These two configurations match perfectly with each other in the 
sense that the 3 magnetic flux quanta of the M~ = 12, i.e., the 3 quasiholes in the center 
of the lower QD "recombine" with the 3 quasiparticles of the = 3 configuration of the 
upper QD. Of course, one cannot label the electrons of different dots due to the inherent 
particle indistinguisibility so one must think in terms of a linear combination of the above 
total configuration with the reversed one: quasiholes in the upper dot and quasiparticles in 
the lower one. Thus, the origin of the u^^^^ = 1 from the spatial inter-dot correlations 
becomes clear. 
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V. CONCLUSIONS 



In this work we have analyzed the new correlated electronic states appearing in double 
quantum dot systems. The simplest case of two electrons was solved analytically By means 
of exact diagonalization of the full interacting Hamiltonian for up to six electrons, minimum- 
isospin ground states were found to appear for short distances between identical quantum 
dots, and to disappear as the distance increased. Mean field calculations revealed the critical 
distances at which the inter-dot correlations were significant, signaling the appearance of 
the new minimum-isospin ground states. 
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FIGURES 

FIG. 1. A schematic picture of a DQDS. Such type of semiconductor structures can be now 
routinely fabricated by means of combined growing and etching techniques. 

FIG. 2. Evolution of the energy spectrum as a function of the distance between QD's of a 
2-electron DQDS with total angular momentum M = 6. Notice the splitting of the / = 1 set of 
states into = and = ±1 states as the distance between QD's increases. The corresponding 
center of mass and relative angular momenta have been explicitly stated in the figure. 

FIG. 3. The phase diagram (angular momentum and isospin) of the 2-electron DQDS ground 
states. 

FIG. 4. (a)Evolution of the absolute GS angular momentum M as a function of the sin- 
gle-particle kinetic energy (magnetic field/confinement) for 4 electrons in the ZDL. (b) The same 
for 6 electrons. Along with M, the values of the total isospin I of the stable states discussed in 
the text are also shown. 

FIG. 5. The excitation spectrum of 6 electrons in a DQDS in the ZDL. The GS is marked 
with a star. Isospin-reversed quasiparticle excitations are distinguished by solid dots and quasihole 
excitations by open dots. The arrows show those subspaces whose GS's will become the absolute 
GS in the presence of isospin (spin) Zeeman energy for increasing B (to the right of M = 15) and 
for decreasing B (to the left of M = 15). 

FIG. 6. (a)Evolution of the absolute GS angular momentum M as a function of the sin- 
gle-particle kinetic energy (magnetic field/confinement) for 6 electrons for D = 50A. (b) The 
same for D = lOOA. 

FIG. 7. Single particle occupation numbers f+ (i^") of the v^^^^ = 1 state as a function of 
the distance. 
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FIG. 8. Evolution of the renormalized GS energy (substracting Neo) as a function of B. Crosses 
show the HF approximation results discussed in the text and dots those from the exact diagonal- 
ization calculations. The inset shows a phase diagram for the inter-dot correlation as a function of 
B and D. 
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